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Abstract 

In this paper we generalize slave-boson mean-field theory for t — J model to 
the time-dependent regime, and derive transport equations for t — J model, 
both in the normal and superconducting states. By eliminating the boson and 
constraint fields exactly in the equations of motion we obtain a set of trans- 
port equations for fermions which have the same form as Landau transport 
equations for normal Fermi liquid and Fermi liquid superconductor, respec- 
tively with all Landau parameters explicity given. Our theory can be viewed 
as a refined version of U(l) Gauge theory where all lattice effects are retained 
and strong correlation effects are reflected as strong Fermi-liquid interactions 
in the transport equation. Some experimental consequences are discussed. 

PACS Numbers: 71.10.Fd, 71.27.+a, 74.20.Mn, 74.72.-h 
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I. INTRODUCTION 



With accumulated experimental evidences it is now generally believed that high-T c super- 
conductors are BCS-like d-wave superconductors close to Mott insulator. It is also believed 
that strong electron correlation does not destroy Fermi-liquid behaviour at low temperature, 
and the ground state of the superconducting states remains Fermi-liquid- (superconductor) 
[1-3] like. Theoretically strong electron correlation is often handled by introducing slave- 
particles, together with a constraint that the sum of electron plus slave particle is equal 
to one on every lattice site in the system. The method allows for approximate mean-field 
solution which gives qualitatively reasonable results in many cases [4,5]. 

Within the slave-particle approach a simple theory that describes a strongly-correlated 
d-wave superconductor is the t — J model in the slave-boson formulation [6]. The slave- 
boson mean-field theory (SBMFT) of the model produces a phase diagram which agrees 
qualitatively with experiments [6]. To avoid complicated mathematics, fluctuations upon 
mean-field theory are often studied in the continuum limit in the language of gauge theories 
[7-9]. Whereas the low temperature phases of gauge theory approaches are Fermi-liquid 
like, non-Fermi liquid behaviour arises rather naturally in the formulation at high temper- 
ature, and can explain qualitatively a number of properties of cuprates [7,8]. However, the 
gauge theory approach disagrees with experiments in some important details [9,10] and it is 
generally believed that the theory is not a realistic representation of high-T c cuprates in its 
present form. 

In this paper we perform a tour de force analysis of Gaussian fluctuations above SBMFT 
of t — J model without introducing continuum approximation. We extend SBMFT for the 
t — J model to the time-dependent regime and study Gaussian fluctuations in the model in 
form of transport equations for both normal and superconducting states. We shall restrict 
ourselves to the zero-temperature, bose-condensed phase in this paper where the mean field 
state is Fermi-liquid like [7-9]. By eliminating the boson and constraint fields exactly in 
the equations of motion we obtain a set of transport equations for fermions of the same 
form as transport equations for normal Fermi liquids and Fermi liquid superconductors 
[11,12], respectively. Effective g-dependent Landau interaction can be extracted directly 
from our transport equations. When computing physical response functions our approach 
gives more refined results compared with U(l) gauge theory [7-9] and offers a precise Fermi 
liquid interpretation to the results. Moreover, with the Fermi liquid representation our 
theory can be extended beyond Gaussian theory more readily. We shall show how strong 
Fermi liquid interactions which become singular in the Mott Insulator limit x — > (x— hole 
concentration) arise naturally as a result of strong correlation in the model and shall discuss 
some of their physical consequences. Unfortunately we discover that the major drawback of 
gauge theory is not remedied by this more complete form of Gaussian theory. 

We begin in section II with a brief revision of slave-boson mean-field theory for the 
t — J model. In section III we generalize SBMFT to the time- dependent regime and derive 
transport equations in the linear response regime. By eliminating the boson (hole) and 
constraint fields exactly we obtain a set of transport equations in terms of fermions only. 
The equations have the same form as corresponding transport equations for Fermi liquids 
in the q — > limit and we can extract (f-dependent Landau interactions from the transport 
equation explicitly. 
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The solutions of the transport equations are discussed in section IV. With a simplified 
form of Landau interaction we compute the density and current response functions and 
show that our refinements does not modify qualitatively the results obtained from U(l) 
gauge theory (See also Appendix A for a direct derivation of gauge theory results). The 
important effects of screening and current renormalization are discussed. The origin of the 
small current carried by quasi-particles j p ~ x [9,13] in gauge theory is interpreted in Landau 
Fermi liquid theory language. Our results are summarized in section V where we comment 
on the lessons we learn and limitations of the theory. 



II. SLAVE-BOSON MEAN-FIELD THEORY FOR T- J MODEL 

We consider a generalized t — J model on a square lattice [6] that includes Coulomb 
interaction between charges [14]. In slave-boson representation the Hamiltonian is 



H = -t J2 (b i b+cic ja + H.C.)-piJ2cic ia + J ]T sJj + ^V^-^Mhbp, (1) 

<i,j>cr icy <i,j> i,j 

where cf a (cicr) and bf(bi) are the fermion (spin) and boson (hole) creation (annihilation) 
operators at site i, respectively. The physical electron operator is f ia = c ia bf. Si = ■ipfaipi, 

where a are Pauli matrices and ipi — ( C *M • The first term in the Hamiltonian represents 



electron hopping where < i,j > denotes nearest neighbor pairs on a square lattice, fi is 
the chemical potential that fix the average density of electrons. The third and forth terms 
represent Heisenberg exchange interaction between electron spins and repulsive Coulomb 
interaction between charges, respectively. The constraint that the number of electrons + 
number of empty site = 1 on every site is introduced in the last term via the Lagrange 
multiplier field \ in the Hamiltonian. 

In SBMFT the various terms in the Hamiltonian are decoupled as follows [6]: 

bib+c+c ja ->< bi >< b+ > ct a c ja + 6*6+ < c+c i(J > - < >< &+ >< c+c i(T >, (2a) 

where < bi >=< bf >=b = ^fx at equilibrium, where x = concentration of holes. 

~l [< AJ > Ajj + A± < A i:j >-< A± >< Ay > (2b) 

~l~ *^ Xij > Xij "I" Xij < Xij > — < Xij X Xij > 5 



Si-Sj 



U 



where < .. > denotes expectation values. = c^Cji — c^c^ and Xij — Ha c ta c jv The A 
term represents decoupling of the Heisenberg interaction in Cooper channel and is responsible 
for superconductivity, whereas the Xij terms represent decoupling in "Fock" channel, and 
gives rise to self-energy corrections to the fermions. The boson (density-density) interaction 
term is decoupled as 

b+bib+bj ^< b > 2 (bf + bi)(b+ + bj), (2c) 
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which is the usual Bogoliubov decoupling [15]. The Lagrange multiplier field is replaced by 
a constant number at equilibrium, 



\ ^< \i >= A. (2d) 

The resulting mean-field Hamiltonian describes two separate species of particles. The 
mean-field Hamiltonian for fermions (spins) describes a BCS-superconductor with 

H f MF = E^t c fc + E [A*(£)(%c_ fel - % c^) + H.C] , (3) 

krj k 

where = -(t6 2 + ^x)t(^) + A - X =< Xij > where j = i + //(// = x,y), and 
7(F) = 2(cos(k x ) + cos(ky)). A(k) = ^-A(cos(k x ) — cos(k y )) represents a fermion pairing 
field with d x 2_ y 2 symmetry, where A =< A i)i+X >. ^fxi(k) is the "Fock" self-energy from 

— — * 

the Heisenberg interaction, whereas A(k) comes from the "Cooper" channel decoupling. The 
mean-field dispersion for the fermion quasi-particles is Ef(k) = + \A(k)\ 2 [6]. Notice 

that the effective hopping from t-term is renormalized by a factor b 2 ~ x in mean-field theory 
and the fermion dispersion is dominated by the "Fock" self=energy. 
The boson (hole) mean-field Hamiltonian is 

H h MF = E «Wp9+ ^ E V®(bt + b-^bf+btf), (4) 
q q 

where V(q) is the Fourier transform of V(r) and e(q) = —txi(q) + A. The bosons are bose- 
condensed in the ground state implying that A = tx7(0, 0). The mean-field dispersion for the 



density excitation is E h (q) = ye(q) 2 + 2xe(q)V(q). For short-ranged repulsive interaction, 
Eh(o) ~ V^q & t small q, whereas for (2D) Coulomb interaction, E h (q) ~ ^fxq (acoustic 
plasmon). The mean-field Hamiltonians together describe a normal Fermi liquid state when 
A = 0, < b >^ 0, x 7^ 0, and d-wave superconducting state when A ^ 0, < b >^ 0,x ^ 
[6]. 

We note that we have followed the so-called U(l) formulation of the t — J model in 
this paper. An alternative formulation is the SU(2) formulation [8]. Mathematically the 
main difference between the two approaches is that the constraint of no-double occupancy is 
written in an SU(2) symmetric way in SU(2) theory, which introduces additional constraints 
Cia c i-a = and its Hermitian conjugate [8]. These additional constraints are automatically 
satisfied if the constraint bfbi + J2 a c ta c iv = 1 is imposed exactly. In mean- field theories 
where constraints are only satisfied only on average, the introduction of SU(2) constraints 
results in more complicated transport equations which are probably more reliable. In the 
low-temperature superconducting states the SU(2) formulation predicts new collective modes 
[16] near momentum q ~ (jr, ir) that are not treated correctly in U(l) theory. However, the 
two approaches give qualitatively similiar results at q ~ (0,0) [16], which is the regime 
we shall concentrate on in this paper. We follow the U(l) formulation in this paper for 
simplicity. 
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III. TIME-DEPENDENT MEAN-FIELD THEORY AND TRANSPORT 
EQUATIONS FOR T - J MODEL 



We now extend SBMFT to the time-dependent regime. The approach can be understood 
by considering the Heisenberg equations of motion (at imaginary time) §-6 = [H,6] for 
operators 2 's that are quadratic in the fermion creation/annihilation operators. Because 
of the presence of four-operator terms in the Hamiltonian, the equations of motion of O2 
will generate terms O4 which involve four-operator terms. The idea of (time-dependent) 
mean-field theory is to decouple the expectation value of the four-operator terms in terms of 
products of two-operator teams in the equation of motion, i.e. < 4 >~< X dtp >, 
thus arriving at self-consistent equations for expectation values of two-operator terms < 
> [11]. The decoupling scheme for SBMFT at equilibrium states will be used in the 
Heisenberg equation of motion. 

We first consider equation of motion for the operator Xij defined in last section in the 
presence of external vector EM field A??*(t) and scalar field ^f £t (r). Following the above 
mentioned procedures, we obtain after some algebra the mean-field equation of motion, 

3 J 

TT- < Xij(r) >=—£" E (< A 5'( r ) >< A fA T ) > + < Xif(r) >< Xj'jir) >) ( 5a ) 
OT 8 j'=i+fi 

3 J 

+Y £ (< A ^'( r ) >< A ^'( r ) > + < ^-'( r ) >< ^'( r ) >) 

i'=j+p, 

-t £ < V(r) >< bf(r) X ^(r) > 

j'=i+fi 

+t ^' t(T) < bjir) >< 6+(r) >< > 

i'=j+fi 

+ (< Xi(r) >-< Aj(r) >) < X ij(r) >, 

where various terms have same meaning as in last section except the introduction of time- 
dependence. The vector field A??*(r) couples to the electrons through the hopping-term [7]. 
Similarly we obtain for the A^- and 6j operators, 

d 3J 

— < A^(r) > = -— J2 (< A *i'( r ) >< Xj'iiT) >-< Xj'i(r) >< A^-(r) >) (5b) 

j'=i+H 

3 J 

+ Y E (< A i*'( r ) >< xtAr) > + < X i'j(r) >< A w (r) >) 

i'=j+jx 

+t E ^ (T) < h(r) X b+(r) X A,,,(r) > 

j'=i+fi 

+t £ e^' (T) < bj{r) >< b+(r) >< A«,(t) > 

i'=j+jx 

-(< Ai(r) > + < Aj-(r) > -2/x) < A -(r) >, 



and 
9 



< &,(r) >= t £ e ^ tW < b Ar) >< Xj>i(r) > -(< A,(r) > +<&r t (r)) < ^) >, 



(5c) 
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where we have coupled the external scalar potential to the bosons [7]. Similar equations of 
motion are also obtained for Aj(r) and bf(r) fields. Notice that there is no equation of 
motion for the < \(r) > field. Instead, it should be chosen such that the average constraint 

< bf(r) >< bi(r) > +J2a < c ta( T ) c i<r( T ) >— 1 is satisfied at all sites i's at all time r. 

To derive transport equations we go to wave-vector space and linearize the equation 
of motion. First we introduce our notations for the Fourier transformed fields. We write 

< Xijij) >— Xij + &Xij{ T )-> where Xij is the equilibrium value of < Xij > an d ^Xij( T ) is the 
time-dependent fluctuation. The Fourier transform field is defined by 

FTiXn] = E e-^'We**-*™* ( Xij + S Xij (r)) = 5 D {t)n % + p^ r), (6a) 

where D = 2 = dimension, is the Fourier transform of Xij an d Pfc(<f> r ) i s the Fourier 
transform of 5 X ij(T~). Notice that (i,j)'s are not restricted to nearest neighbor sites in the 
Fourier transforms. It can be shown from the definition of X ij that p^(o) = P%{—q)- Similarly 
we define for the A^- field, 

FT[Aij] = E e-iik+m-n^k-m-r, + 5A . . (r) ) = 5 D (q)A^ + A R (q, r), (6b) 

where A^ is the Fourier transform of A^. It can be shown that A^(q) = A_g(q). For the 
boson and constraint fields, we define 

FT [H^)i\ =E eif?i (b(X) + 5bi(T)(5Xi(T))) = 5 D (q)^(X) +feq-(r)(A_q-(r)). (6c) 



Next we linear the equations of motion (5). We consider small time-dependent fluctu- 
ations 5 X ij{T~), 5Aij(r), etc. and keep to linear order time-dependent fluctuations in both 
sides of equations (5). Introducing Fourier transforms via equation (6), we obtain after 
some lengthy but straightforward algebra, six linearized equations of motion for the fluctu- 
ating fields p%{q,r),p+ n (-q,T), A s (q,r), A+-(-£ t), &_ ? -(r) and fc±(r). The A g -(r) field will 
be determined by the constraint equation. We shall simplify our notation by not displacing 
the explicit time-dependence of the fields entering the transport equations in the following. 
First we consider the equations of motion for the fermion fields. 

The equations of motion for the p%(q) and Ag(g*) fields can be written in a matrix form. 
Introducing the vector notation 

( pM \ 

VA^(-g)/ 



we obtain after Fourier transforming in time, 



O J 

G^qM^M = — J J2 cos(A; M ) (AM^+hm;(-Q)) (?) 

fi=x,y 



X r7 - £cos(* M ) (2fcos(|)5(&_*+&f) + ^(fl 



+HI(q) Esm(^) (2t sin(|)6(6_ f - b±) + 2ifb 2 Af(q) - ^(q)] 
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where 



I iuj 



V 



^k+q/2 + £j 





■k-q/2 



MA, 
8 ^k+q/2 

-MA- 

8 ^k-q/2 







^k-q/2 + ^k+q/2 

MA. 
8 ^k-q/2 

8 ^k+q/2 



MA. 

8 ^k+q/2 

MA. 
^k-q/2 

+ £fc-g/2 







-MA- 

8 ^k-q/2 

_MA- 

8 ^k+q/2 





= 



'^k+q/2 
'^k-q/2 

i 





n fe+g/2 + n fc-g/2 



Aj(9) = 



A fc+<?/2 





and 



/ U k-q/2 
k 

-A 



fc+g/2 

n k+q/2 ~ n k-q/2 



k+q/2 ^k-q/2 
V ^k+q/2 + ^k-q/2 ) 



( n k- 


-<f/2 


- n fe+g/2 


\ 


n k- 


-9/2 


_ n fe+<f/2 






-9/2 


~ ^k+q/2 




\A fc - 


-9/2 


~ ^k+q/2 


J 



where p^(q) = Eg cos(/c M )p^(g*), = sin^p^g) and A M (g) = E^ cos(/c M )A^(g). 

Writing A iji+ ^ = s^A^i + |)|e"^ + ^, where s^y) = +(— )1, we obtain with Eq. (6b) for 
small fluctuations A iji+ £ ~ A^ + s M 5|A^|(i + f ) + iA^^i + f ), and 



A„(g) =s„<J|A M |(g)+iA^(g), 



(9a) 



showing that the real and imaginary parts of A^(q) represent the amplitude and phase 
fluctuations of the superconducting order parameter, respectively. Similarly, writing Xi,i+ft — 
IxAii + f )e ia ^ t+ ^\ we obtain for small fluctuations, 



(9b) 



showing that p s ^ a '(q) represents the amplitude (phase) fluctuations of x- Notice that a^q) 
is identified as an internal gauge field in gauge theory [7]. It is also proportional to the 
mean-field fermion paramagnetic current, which is given by jZ(q) = + 2t6 2 )p^(g). 

The structure of the fermion matrix transport equation is quite similar to transport 
equations for usual Fermi liquid superconductors if only spin-singlet excitations are con- 
sidered [17]. The absence of spin-triplet excitation in the matrix transport equation is a 
consequence of SBMFT which consider only decouplings in spin-singlet channel. is 
the standard kinetic term for Fermi-liquid superconductors [17] and the coupling to A and 
ps(a) figi^g j n the equations of motion can be understood easily from standard many-body 
theory for fermions. They are required in a conserving approximation in treating fermion- 
fermion interaction [18]. In equilibrium SBMFT, the Heisenberg interaction is decoupled 
in two channels resulting in the Cooper pairing field and "Fock" self-energy. A conserving 
approximation requires that corresponding vextex corrections must exist when computing 
the particle-hole and particle-particle propagators. They are precisely the A and p s(a ) terms 
in the equations of motion. For ordinary Fermi liquid superconductors inclusion of supercon- 
ductor order parameter fluctuations is essential to restore the correct longitudinal density 
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response for superconductors [19] and inclusion of vertex correction coming from self-energy 
is essential to obtain the correct Fermi liquid interactions [18]. Similar results are obtained 
in SBMFT. The main difference between SBMFT and ordinary Fermi liquid superconductor 
is the presence of boson components an constraints in the system so that in addition to the 
above vertex corrections, fermions in SBMFT also couple to the boson and constraint fields 
in the equations of motion. 

The two linearized equations of motion for the boson fields are 



H=x,y 

and 



V V 

v 

+2itbx E M^)Af{q)-K\^+^ ex \-q)) 

H=x,y 



fb 

{-icv - e(q))bt ? = PV(q)(b,+ b + _,) - - E7(fc' + ?/2)^(5) (10b) 



V V 

k> 

-2ztbx E sm(^Mf($)+6(A_f +$«*(-$)). 

H=x,y 



It is convenient to introduce the density-phase representation for bosons 6, = riie l6i . 
Using Eq. (10) we obtain for small fluctuations, b(b-$ + bt) = n(—q) and b(b_^ — bt) = 
2fb 2 9(-q), 

n(-q) = J-^B^\ f 2 *£ cos (f " (A*+ (Ha) 



and 



(iou) 2 


-Eh® 2 ) \ 




2b 2 {iuj) 


1 fa\ 


) 2 -E h (q) 2 ) 




-2(iu) , 



2<«(-5) = ( ^ )2 _^ 2) (2tE-s(|) P ^(g1 - (A, + $-(g))j (lib) 

In the linear response regime, the constraint field A^ should be choosen such that the 
density fluctuations of the fermion field is exactly balanced by the density fluctuations of 
the boson field, i.e. p(q) = J2kPk(o) = —n(—q). With eq. (11a), we obtain for the A field, 

A ? - = -<& ext (q) + 2t E cos (y K(<?) " VT^—A^) (I 2 ) 



-^)2tEsm(|) {pt(q)+ixA?(q)), 



where Xh(q, iw) = ui^^Jfff ls ^ ne mean-field density- density response function of the holes 
[7]. Putting eq. (12) in Eq. (lib), we obtain 
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2tbh(q)6(-Q) = -4^E«in(f ) (ft® +ixA?(qj) , (13) 

which is the continuity equation. To see that we rewrite Eq.(13) as 

^p(g)+X:2^in(|)J M (g) = 0, (14a) 



where 



JM = 2bHx (2 sm(^)9(-q) - a,(q) - Af{q)) , (14b) 



is the physical current. The current expression (14b) can also be obtained directly from the 
current expression derived from t-term, where Jij ~ —t J2 a {^ Al] c ta c j^i^ ~~ c - c -) • m mean- 
field approximation, ~ —t £ CT (j> 2 (< cf a c j<T > -c.c.) + ^fo 2 (e l(e<_ ^ +j4 ^* ) - c.c.)), which for 
small fluctuations is precisely Eq. (14b) noting that < cl a Cj a >~ xe iaij . 
Notice that using Eq. (13), we can also write 

= -$ ext {q) + 2t E 008(^(5) + §|| + iMe(-g). 

Together with the constraint equation p(q) = —n(—q), we can eliminate the boson density 
and constraint fields from the fermion equations of motion. The coupling to boson phase 
9(—q) can be eliminated through a gauge transformation, c ia — > c iu e~ l6i and bf — > bfe l9i 
which absorb the boson phase on the fermions. It is straightforward to show that under this 
transformation, the fermion field ^^{q) in the transport equation transform in the linear 
response regime as follows: 

PM PM + ifak-w - n k+q/2) 9 (-V)i ( 15 ) 
A^q) -> -i(A s _y 2 + A Uql2 )0(-q), 

a^q) -> a^q) - 2sin(-^)0(-g), 

^(0^(0" 2 cos(f )*("$)■ 

Together with Eqs. (15) and (7), it can be shown that the gauge transformation removes 
all the 9(— q) fields in the equations of motion for the transformed fermion fields. Putting 
together we obtain, 

G£{i!M*M = -¥E«»(*/.) {hmM + A^1A;(-g)) (16) 



4 , 



+HKq) 



+ X>s(^) (2tcos(|)p(g) - ^(g) 
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where all terms have the same meaning as before except they are now expressed in terms of 
the transformed fermion fields. Notice that Eh(q) 2 / (2b 2 e(q)) = e(q)/2b 2 + V(q). 

Equation (16) is the main result obtained in this paper. We note that after eliminating 
the bosons the transport equations for the fermions have the same form as a (<f- dependent) 
transport equations for ordinary Fermi liquid superconductors. The dynamics of the bosons 
is completely eliminated reflecting their 'slave" nature. Besides coupling of fermions to 
superconducting order parameter fluctuations, the fermions couple to each other through 
Landau molecular-field type interactions. We can extract effective g-dependent Landau 
interactions by comparing Eq. (16) with Landau transport equations for superconductors 
[12]. We obtain, 

h f (q) = ^ (§r + v (q) - x E cos (^ - p») + 2* E cos(|)(cos(Ag + CO s( Pm ))) , (17) 

which corresponds to a rather regular Fermi liquid. Notice that the density-density Landau 
interaction diverges in the b 2 — > limit, which we shall see indicates the vanishing of 
compressibility in the insulating state. The effects of strong correlation will be discussed in 
more detail in the next section. 

We note that with Eq (17) we can also introduce quasi-particle scattering into the trans- 
port equation through usual Fermi liquid phenomenology [11,15]. The life-time of a quasi- 
particle with momentum p is given by the Boltzmann equation expression [11], 

- ~ 2vr £ £ W(p, k; q)5(E f (p) + E f (k) - E f (p + q) - E f (k - q))n~ k {l - n^)(l - n^), 

T P q k 

— * 

where the scattering rate W(p, k; q) can be identified as the absolute value square of the 
scattering amplitude \A(p, k;q)\ 2 and can be computed with our Landau interactions (17) 
[11]. We recall that our transport equation is derived only for spin-singlet channel, and 
spin-triplet landau interactions are missing. 

For ease of analysis it is convenient to separate the coupling of fermions to 'longitudinal' 
and 'transverse' fields explicitly. To achieve that we work in the Coulomb gauge such that 
A^^q) is purely transverse and we write 

"M =<(?)+ (18a) 

where p^(q) = ix a ti(o) an d 

We then perform a gauge transformation to rewrite the coupling of the fermion fields to 
a l ^(q) in terms of a coupling to scalar field A|. Using the continuity equation and keeping 
in mind that the boson phases are already absorbed in the fermions in Eq. (16) we obtain 
J2n 2? sin(-j^) Jfj,(q) = — 4db 2 tx J2u sin(^ t )a /i (g) = icup(q). With the gauge transformation 
relation 2sin(^)A°. = (1 — z)i(iu)a l fl (q), we obtain = —(1 — z) ^jr^ p(o), where 



^ sin(^ 



sin(f; 



E,[sin(f)] 2 

sin(f)sin(f) ' 
E,[sin(f )] 2 , 




(18b) 



a v (q). 
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2tb 2 

Putting this back into the Equation of motion (16), we obtain an matrix equation of 
motion for fermions where the longitudinal and transverse interactions are separated, 



(19a) 



-^V,1 + 2.Ecos(|K(,1-^^^ 

3J 



+ £cos(*„) (2tcos(|)p(g) - -p^(q) 



where 



2b 2 e{q) 



(l-z)(iu) 2 -E h (q) 2 ' 
is a renormalized density-density response function for the holes. 



(19b) 



IV. SOLUTIONS TO THE TRANSPORT EQUATIONS 

In this section we shall study the solutions of the transport equations in both the normal 
and superconducting states, concentrated in the small wavevector q limit and neglecting the 
quasi-particle life-time effects. We shall first consider the density and (transverse) current 
responses in the normal state and shall compare our results with corresponding results in 
U(l) gauge theory [7]. We shall then study the solutions of the transport equation in more 
details where quasi-particle properties will be examined. A similar study will then be carried 
out for the superconducting state. 

The full equations of motion for fermions involving the complete Landau interactions 
(17) are too complicated and cannot be solved exactly. To have a qualitative understanding 
of the solutions we approximate the Landau interactions by keeping only density-density 
and current- current interactions, i.e. we approximate 

~ \ {U £ff {q) "jE sin (M **(pSj > (2°) 

where U e ^(q) ~ ^2 V(,0) d(q)t — bj , a, and b constants of order 0(1) obtained from 
averaging the third and fourth terms in the Landau interactions (17) over the fermi surfaces. 
The effect of the at—bJ term is to renormalize the charge-charge interaction V(q) by adding a 
repulsive short-ranged potential of order ~ t. After making this approximation the transport 
equations become 
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G£(£ = -^E «»(*/«) (A^^A^g) + A^A^-g)) (21) 



4 „ 

- — L , p(g)) 

where x RR (Qi ^ U} ) is given by Eq. (19b) with a renormalized charge-charge interaction 
V R (q) ~ V(q) + at — bJ . We shall study this simplified matrix transport equation in 
the following. 



A. solution in the normal state 

In the normal state the fermion matrix equations of motion become a single equation 



® ext (-q) - mk ■ \ P(o) ( 22 ) 
x RR {(iM 



f 3 T 

+ £sin(A;,) hitb 2 A^{q) - —ixa^q) 

In the small <f limit the longitudinal and transverse responses are decoupled from each 
other and the density and current response functions can be solved easily. We obtain for 
the density-density response function, 

x(, '" )= *» + rt»' (23) 

where x d (g, w) is the density-density response function of the system and xf(g, cj) is the 
corresponding free fermion response function, i.e., xtiQi 00 ) = J2kXok(k,u), where 



To compute the current-current response function, we note that the physical current 
(14b) can be written as J^q) = —zjj^ = —z(j^ + jt) after the bosons are eliminated, where 

3u = (t 2 ^ 2 )x°m(^) an d if — (2t6 2 + ^xlx-^^io) are the fermion paramagnetic and 
diamagnetic currents, respectively. With this we obtain from solving equation (22) 

x [q,u — 7; \ tt=. — ^— — rr^ — v 

{l-z)x t s {q,uj) + x t h {q,uj) 

where x*(?> ^) is transverse current-current response functions of the system, and xKQi 00 ) = 
+ 2t6 2 ) sin(/c /i )] 2 x ^(?, is the corresponding free fermion response function. 
Xh(q,uj) = —2b 2 tx- We have written our results in a way which can be compared with 
gauge theory results easily [7,8]. It is obvious that the transport equation produces essen- 
tially the same result as Ioffe-Larkin result for the density and current response functions, 
except the additional z factors and the renormalization of density-density interaction V R (q). 
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We now examine the solutions of the transport equation in more detail. Following Fermi 
liquid theory, we consider solutions of the transports equation as eigen-excitations of the 
system [11,20]. There are in general two kinds of solutions, dressed particle- hole excitations 
and collective modes [11]. With our simplified Landau interaction, there exists only collective 
modes associated with density fluctuations in the system which can be determined from the 
poles of the density response functions. At small q, the dispersion of the collective mode goes 



as ujq- ~ y2b 2 e(q)(V R (q) + A^(0) _1 )/(1 — z), where iV(0) is the density of states on Fermi 
surface. Writing V(q) ~ tv (q) where v (q) is dimensionless, we obtain for small x oo^ goes as 
ty/qx for Coulomb interaction, and goes as tq\fx for short ranged interaction. Well defined 
collective mode exists at \q\ < x(t/J) 2 for long-ranged Coulomb interaction. Otherwise the 
collective mode is overdamped by the particle-hole continuum. 

We next consider solutions corresponding to dressed particle-hole pair excitations. In 
the normal state, we look for solutions of form [11] 

pM = pf(q) + s Pk(q), 

where P^\q) = S(uj — ^+q-/2 + £,k-q/2)d D (P ~~ corresponding to a bare particle-hole pair 
with relative momentum p and energy Cp+q/2 — Ck-g/z • Putting this into the equation of 
motion (22) we obtain 



pn(<t) = -pf (o) + x M UJ ) 



1 / 3 7 

® eXt (~fi + vflflr*,-, ^ (g) - £sin(^) (2itb 2 Af(q) - -ixal{q) 

(25) 



The properties of the quasi-particles can be determined from solution of Eq. (25) if we 
identify p%{f) as a local quasi-particle density as in Fermi-liquid theory, it is easy to solve 
Eq. (25) to obtain 



Xh R {Q, - Z,k-q/2 



= RR( ^ : : at — * — -/ w> w 

Xh KQ, ?P+<f/2 - Sk-q/2' +Xs{q, Kp+q/2 ~ Kk-q/2> 

where p(q) = X^Pfc(<f) ls t ne total charge carried by the dressed particle-hole pair [11] and 
p(°) (g*) ~ 1 is the charge carried by the bare particle-hole pair. Overall speaking we find that 
the charge carried by the eigen-particle-hole pair (and the corresponding quasi-particle) is 
renormalized by a factor of order \ RR ~ x from its bare value because of screening [11]. For 
long-range Coulomb interaction the charge carried by the dressed particle-hole pair vanishes 
in the small q and u = limit. Similarly, we can also compute the (transverse) current 
carried by quasi-particle from Eq. (25). We obtain 

m = tt^ w-l ^- v J m (^ (27) 

where J ( -°- )t (g*) is the transverse current carried by the bare particle- hole pair. As in charge 
response, the overall current carried by the quasi-particle is renormalized by a factor of 
order Xh ~ x - However, in the small q and uo = limit, the free fermion transverse current 
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response function xl vanishes because of gauge invariance whereas Xh remains finite ~ x, 
and the transverse current carried by quasi-particle is unrenormalized from its bare value. 

Our results can be understood readily from a Fermi liquid picture. With our simplified 
Landau interaction the long wavelength properties of the system are governed by the effec- 
tive mass m*/m and two Landau parameter F ~ V R (0)N(0) (for short-ranged interaction) 
and Fi = (-2—1). The effective mass correction is equal to ^ — fb2+ ^ J ^ g - ) - , which is of order 

t/J. Notice that the Fermi liquid identity 1 + Fi — — which is rigorously satisfied for trans- 
lationally invariant systems [11,12] is not satisfied here. Microscopically F\ is coming from 
the vertex correction corresponding to the "Fock" decoupling of the Heisenberg interaction 
and vertex correction associated with slave-boson renormalization of bandwidth t — > b 2 t is 
absent in SBMFT. 

It is straightforward to see that our density-density response function is consistent with 
Fermi liquid theory. The static density-density response function is given by 

X-(?,0) = : 



u°ff(q)xi(q,oy 



in agreement with Fermi liquid theory [11]. Notice that x d (<f)0) goes to zero in the Mott 
insulator limit x — > 0, (lq) 2 » 1, where / = 1/y/x is the average distance between holes, re- 
flecting the vanishing of compressibility in the insulator state. Notice that external impurity 
potentials V A/ (g*)'s are screened by the factor 1 — U e ^{q)x'l 0) correspondingly, implying 
that the effect of impurity scattering on quasi-particles is weakened by strong correlation 
and vanishes as the Mott insulator state is approached [14]. 

For the current response, the current expression (14b) implies that the physical current 
operator written in terms of the fermion current operator is given by J^q) = (1 + Fi)j^(q), 
in agreement with Fermi liquid theory [11,12]. Our results for the current response can be 
obtained from Fermi liquid theory with this identification [11,12]. Notice that the "back- 
flow" current represented by F\ exactly cancels the "bare" current in the Mott-insulator 
limit (z ~ x — > 0), implying that the quasi-particles carry zero current in this limit, which 
seems to be consistent with an insulator phenomenology. 



B. solution in the superconducting state 

Next we analyse the fermion matrix transport equation in the superconducting state 
with the approximate Landau interactions (20). 

As in the normal state we study electromagnetic responses of the system in the small q 
limit where the density and transverse current responses are decoupled. First we consider 
the density response. With the simplified form of Landau interaction, the matrix transport 
equation becomes 

O 7 

g-^m^m = - t e««w (hmM + hmu-v)) ( 28 ) 

X KK (q,iu) 
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which represents responses of a "pure" d-wave superconductor with no Landau interaction, 
to an effective external potential — <J> e:ct (— q) — ^m^jjj Pio)- Defining Xds(& u ) as the density- 
density response function of the pure d-wave superconductor, we find that the density- density 
response function of the system u;) is again give by equation (23), except that xf(<f, u;) 
is replaced by the corresponding response function xls(<L u ) °f a d-wave superconductor. 
Similar conclusion is also reached for the transverse current response, with xKQi 00 ) replaced 
by the superconducting counter part u ) 1X1 Eq- (24). These results are consistent with 

a generalized U(l) gauge theory for the superconducting state [7,8]. 

As in U(l) theory, the superfluid density deduced from x*(0,0) is given by p s ~ 
a'x — b'x 2 T at small x and low temperature T, where a' and b' are numerical factors of 
order 0(1) [8-10,13], and the b'x 2 T term is in disagreement with experimental results on 
the London penetration depth. The additional (1 — z) factor that appears in our transport 
equation solution does not modify this conclusion qualitatively. The x-dependence in the su- 
perfluid density originates from the Fermi liquid renormalization of the quasiparticle current 
Jn(q) = (1 + Fi)jf(q) and is a fundamental property of SBMFT [11,12]. The disagreement 
with experiments suggests a fundamental problem associated with SBMFT for the high-T c 
cuprates [13]. 

To study charge response, we neglect the small dissipative effects coming from nodal 
fermions and assume 

d Kg 2 
X*(9' W ) = u,*-Kuq*> 

where K,u ~ J and are the phase stiffness and inverse compressibility for the "pure" d- 
wave superconductor, respectively. We obtain from Eq. (23) that the K and U factors are 
renormalized by Landau interactions at small q to K — > (2xK/(2x + (1 — ~ %t and 

u — > u + V R (0) for short-ranged interaction V R (q), in agreement with U(l) theory [21]. For 
Coulomb interaction the acoustic plasmon mode remains in superconducting state at small 
enough q. Notice that the strong renormalization of compressibility and screening of impurity 
potential n the limit {lq) 2 » 1 remains qualitatively the same in the superconducting 
state and may be the explanation why the nanometer-scale gap inhomogeneity revealed by 
STM experiments on BSCCO surface [22] does not seem to introduce strong quasi-particle 
scattering effect [14]. 

The quasi-particle properties in the superconducting state can be studied from the trans- 
port equations as in the normal state. We look for solutions of form ^^(q) = ^f^\q)+5^^(q), 

where ^f(q) = 5(u-E f (p-q/2)-E f (k-q/2))5 D (p-k)a p (q), where £/(£)'s are mean-field 
quasi-particle energies and 

u p-q/2Vp+q 
—Vp~-q~/2Vp+q 
V Up--ql2Up + q j 

where u 2 - — h(l + and v 2 - — \(1 — are superconductor coherence factors. The 

solution represents a dressed pair-excitation in the superconducting state. Proceeding as 
before we find that the charge and transverse current carried by the dressed particle-hole 



a p (q) = 
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pair are given by Eqs. (23) and (24) as in the normal state, except that the normal state 
fermion response functions xt an d xl are replaced by the corresponding functions in the 
superconducting state Xds an d Xdsi respectively. In particular, we observe that because of 
vanishing of paramagnetic current, the transverse current carried by a quasi-particle is now 
renormalized by a factor z ~ x, which results in the b'x 2 T dependence of the superfluid 
density. The origin of the z factor can be chased back to the current renormalization 
Jn{q) = (1 + Fi)jf(q) as we have discussed. 

Notice that besides carrying charge and current, quasi-particles also carry with them 
disturbances in the mean- field order parameters Xij an d A y -. It is found in SU(2) theory that 
low energy collective modes associated with fluctuations of the Xij andA y - order parameters 
exist at momentum around (tt, it) and may enrich quasi-particle properties considerably [16]. 
These effects are not considered with our simplified form of Landau interaction (20). 

V. CONCLUSION AND COMMENTS 

We now summarize our findings in this paper and discuss a few questions associated with 
our theory. 

By deriving complete transport equations we show in this paper how Gaussian fluctua- 
tions described by SBMFT oft — J model can be understood in the language of Fermi liquid 
theory. An effective Fermi liquid Hamiltonian is derived which form the basis of our calcula- 
tion of response function and quasi-particle properties. Our theory provides a more accurate 
description of Gaussian fluctuations in t — J model compared with U(l) gauge theory and 
has the advantage that calculations beyond Gaussian theory (e.g. various transport coeffi- 
cients) can be conveniently performed within the framework of Fermi liquid phenomenology. 
Unfortunately the main drawback of gauge theories associated with quasi-particle current 
renormalization is not remedied in our treatment. Notice that our formulation can be gener- 
alized to the high temperature phase where bose-condensation is absent in mean-field theory 
with a modified mean-field decoupling scheme where < b { >= but < bfbj >, < bfij >^ 0. 
It will be interesting to find out whether the high temperature phase can be identified as 
some kind of Fermi Liquid phase at high temperature in this approach. 

A more fundamental question is to examine general features in our theory which are not 
specific to the t—J model. To be specific we consider a lattice fermion model with hopping to 
both nearest and next-nearest neighbor sites and with (repulsive) density- density and spin- 
spin interactions different from the t — J model. The only common feature of the model to 
t — J model is the presence of strong-correlation handled by introducing slave-bosons. It 
is not difficult to see that the equations of motion for bosons will be essentially the same 
as the corresponding equations in t — J model, except that the mean-field dispersion e(q) 
and density-density interaction V(q) are different. For repulsive interactions the behaviour 
of bosons will be very similar to bosons in t — J model. The expression for Ag-field in linear 
response regime will also be very similar to those in t — J model, since it is insensitive 
to the exact form of the spin-spin interaction. Therefore the main difference comes in the 
fermion part where a large variety of mean-field phases are possible, depending on details of 
the spin-spin interaction. After elimination of boson phase a generalized Fermi-liquid type 
description for fluctuations around the fermion mean-field state will be obtained for the 
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lattice model, the precise forms of the Landau interactions will depend on the microscopic 
Hamiltonian and the mean-field fermion state. 

Notice however we expect that the Landau interactions will have some common features 
as in t — J model. First of all, we expect that an density-density interaction term of 
form f£-(q) ~ E h (qf /2b 2 e(q) will remain, since it comes from eliminating the bosons and 
constraint fields and does not depend on the fermion field. In fact, the presence of this term 
guarantees that the system becomes incompressible in the x — > limit. Secondly, we expect 
that current-current interaction term Fi ~ ax — 1, where a is a numerical constant of order 
0(1) will remain, as long as the fermion mean-field bandwidth remains finite in the limit 
x — > 0. The existence of this factor in SMBFT can be inferred from mean-field decoupling 
of the hopping terms which always introduce a factor ~ x in the coupling to external EM 
field. As a result quasi-particle carries vanishing current in the Mott Insulator limit x — > 0. 
Notice that these are the two main features that dictates the density and current response 
of the SBMFT treatment in the small q limit. The temperature dependence of London 
penetration depth in high-T c cuprates suggests that current carried by quasi-particles on 
the Fermi surface does not vanish in the x — > limit. It would be interesting to examine 
whether this feature is unique to high-T c cuprates or whether there exists systems where 
quasi-particle current indeed vanished continuously when approaching the Mott-Insulator 
state. 

Finally we note that U(l) formalism is believed to be insufficient when describing prop- 
erties associated slightly broken SU(2) symmetry. For simplicity we have concentrated at 
region q ~ (0, 0) in this paper where SU(2) effect is not strong. To incorporate SU(2) 
fluctuations correctly a SU(2) version of transport equations is required. 
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APPENDIX A: IOFFE-LARKIN FORMULA FROM TRANSPORTS EQUATION 

To make direct contact with U(l) gauge theory [7] we consider the equations of motion 
before elimination of bosons. First we consider the current expression (14b). We note that 
the current expression can be interpreted as a boson (hole) current under the influence of 
effective gauge field a^q) + A e * l (q). Alternatively we can rewrite the current expression as 

O J 

JM = "(2^ + T X)X (aM ~ (1 - + * - , (Al) 

where a^(g) = 2 sin(— q^/2)9(— q) is an effective gauge field coming from boson phase fluc- 
tuations. The current expression can be interpreted as (-ve) a fermion current under the 
influence of effective gauge field (1 — z)a /x (q) — z{A e l ^ t {q) - a b ^(q)). The existence of two 
interpretations for the physical current is a basic feature of U(l) gauge theory. 

We shall consider the normal state for simplicity. The superconducting state can be 
treated in a similar way. In this case = and the fermion matrix transport equation 
is reduced to a single equation for the p%(q) field. We further neglect direct coupling of 
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the bosons and fermions to the p s ^{q) and n(—q) fields in these equations and compute the 
density and current responses in both boson and fermion current interpretations. Using Eq. 
(11) we obtain for the density and transverse current in "boson" -representation (in Coulomb 
gauge), 

n(-q) =Xh(q,u)(\ i ?+\ a f+® ext (q)), (A2) 

where A| comes from gauge transforming the longitudinal part of the a^(q) field. a^(q) is the 
transverse part of a^(q). Using Eq. (7), we obtain for the density and transverse currents in 
" fermion" -representation, 

p($=X*tfM(>*+0- -*)>$) (A3) 
JM = -Xl(q, <w)((l - zK(q) + zAf(q)), 

In the z — > limit, we can neglect the O(z) terms in equation (A3) and the two equations 
can be solved self-consistently easily. We obtain the well known Ioffe-larkin expression for 
the density and current response functions, 



Xh{q,u)+x d s {q,uY 
and 



Xhfav) +x t s (q,v)' 

A parallel analysis can also be done for the superconducting state. 
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